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any deposit (1). When constructing variograms, the samples must belong to one zone 
to represent the spatial variability of that zone (2). These zones can be defined by 
studying the geology of the ore body including its structural, mineralogical, and 
chemical properties. The database is then created for a specific zone, in which all 
segments are assumed to have the same geological features. 

Based on previous geological studies, which were carried out on the Ghorabi iron ore 
deposit, it is clear that the orebody could be divided structurally and stratigraphically 
into four zones. This division might affect the mineralogical characteristics of the 
orebody and the ultimate ore reserve calculation. However, this could be investigated 
through geostatistical modeling of the orebody by considering its geology and hence 
variograms have been constructed for individual different zones. 

For sake of comparison, it is then assumed that the structural geology has no effect on 
the mineralization of the orebody, and the deposit whole area is considered as one 
zone for which variograms could be constructed. 

To clarify the advantage of geostatistical analysis over conventional statistical 
analysis, the distribution of iron content and its parameters; mean, standard deviation 
and coefficient of variation, within these zones has been performed statistically. This 
statistical analysis could refer to different characteristics of iron content within the 
different zones. 

GEOLOGY OF GHORABI OREBODY: 

Gabal Ghorabi is located on the trace of a giant dextral right stepping wrench fault 
system that extends further in the northeast direction. In the mean time the master- 
fault affecting Gabal Ghorabi makes two levels of iron ore, one high to the north** 
occupying the foot-wall while the lower iron ore level is situated in the hanging-wall 
( 3 ). 

Ore geology and ore genesis : 

The orebody could be classified into the following types (3): 

(a) Stratiform Cenomanian ironstones: 

This ore type forms stratiform layers and lenses of variable thicknesses, confined, to 
the upper mudstone dominated member of the Cenomanian Bahariya Formation. The 
ore is of muddy nature and ore layers and lenses often intercalated by glauconitic 
mudstones and siltstones. 

(b) Stratiform to stratabound Middle-Eocene ironstones comprising: 

- Stratiform oolitic or oolitic-nummulitic ironstones (3-6 m in thickness) deposited 
directly on. the Cenomanian elastics and the associated ironstones representing 
formation under shallow marine condition (3). 

- Stratabouned surfacial duricrust developed oolitic-pisolitic ironstones which form 
the hard cap demarcating the summit of gabal Ghorabi. The ore is of brecciated nature 
dominated by crustified textures or cave fillings formed as a result of subaerial 
denudation of pre-existing iron rich sediment. 

STRUCTURAL MODEL: 

Further north, the northeastern plateau is affected by another giant right-lateral fault, 
known as the Ghorabi fault system (Fig 1). This fault system extends through the 
Eocene sequences to the north of the Depression. It brought the Cenomanian elastics 
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samples with respect to each other are not taken into account (2). This could be 
considered a reason for the coefficients of variation not to give any clear difference 
between the studied zones and hence reflects the insufficiency of traditional statistics 
in analyzing the variability of Ghorabi orebody. 

Geostatistical analysis, which allows the relative sample positions to be taken into 
account, car* therefore be used to check this supposition. 

CONSTRUCTING VARIOGRAMS: 

Construction of an experimental variogram is the first step in any geostatistical 
analysis. It can be computed from a set of randomly spaced data through finding pairs 
of data that are oriented in the required direction, determining the distance between 
the samples, then summing the squared differences of the grades and dividing by the 
number of pairs (1). 

Constructed variograms reflect the variability of the ore mineralization and could be 
used to differentiate between characteristics of the different zones of an orebody and 
as a base of the next steps in the geostatistical evaluation, such as kriging. To illustrate 
the effect of geology of the orebody, directional variograms have been constructed for 
each zone and for the whole area to check the directional anisotropy characteristics for 
each case and for the whole area. 

TechBase software, which is a system designed to aid geologists and engineers in 
mine evaluation and mine planning, was used to construct the variograms. According 
to the data arrangement, the required parameters for constructing variograms were 
introduced to the program. A lag interval of 100 m with 10% tolerance was^' 
considered when constructing horizontal variograms. When the directional variograms 
were constructed, the azimuth (bearing, clock-wise from north) was considered as 0° 
for N-S, 90° for E-W, 45° for the diagonal direction NE-SW and 135° for NW-SE 
with a 22.5° tolerance for any direction. 

Directional anisotropy: 

Directional variogram reveals the change in the variogram parameters as the direction 
changes (7). The directional variograms for the whole area as shown in Fig. 4 
demonstrate the absence of anisotropy when Ghorabi iron deposit is considered as one 
zone and the geological characteristics are neglected. 

The variogram curves are sufficiently close to each other through the different 
directions and appear to be the same. Also, the fitted spherical models through the 
four directions, as shown in Fig. 9, do not show any clear differences. This 
information is also shown in Table 2. The orebody of Ghorabi area, at this poiht, 
could be considered isotropic when neglecting the effect of its geologic structure. 

However, when considering the effect of the geology of the orebody, i.e., considering 
both of faults and stratigraphy of the orebody, directional variograms within the 
different zones reflect the presence of directional anisotropy through the different 
-zones. The appearance of linear model on some directional variogram reflects high 
continuity through these directions, i.e., the ore mineralization is considered to have 
low -variability. Direction NE-SW through zones (S) and (E), E-W of zones (E) and 
(C),'N-S'bf zbiie (N) give higher continuity than other directions. Also, the low values 
of C 0 and C on the’ ipherical ’models of the other directional variograms illustrate that 
the ore mineralization varies from one direction to another within narrow range. 




Zonal anisotropy: 
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Table (1): Statistical Parameters 



Zone 


mean 


S.D. 


Cov 


N 


46.97 


2.8 


0.06 


E 


47.07 


3.9 


0.08 


C 


46.42 


4.4 


0.09 


s 


But 


3.1 


0.06 


Ghorabi 




4.00 


0.08 



Table (2): Directional variogram parameters 



ZONE 


N-S 


E-W 


NE-SW 


NW-SE 


Ghorabi Co 


0.5 


0.7 


0 


0 


C 


11.9 


11.4 


14.8 


11.3 


a 


320 


320 


300 


300 


(N) Co 


0 


0 


0 


1.07 


C 


without 


12.5 


7.0 


12.4 


a 


without 


300 


220 


250 


(C) Co 


0 


I.l 


0.6 


0.75 


C 


8.9 


without 


11.8 


14.6 


a 


300 


without 


300 


300 


(E) • Co 


0 


0 


0 


0 


C 


15.6 


without 


without 


11.7 


a 


300 


without 


without 


250 


(S) Co 


0 


wmm 


1.4 


0 


C 


9.9 




without 


9.0. < 


a 


300 




without 


300 



Table (3): Zonal variogram parameters 



ZONE 


Range of 
influence 


Sill 


Nugget 

effect 


Ghorabi (T) 


400 m 




. * 1 


(E) 


350 


B 


0 


(N) 


mm 


7.4 


0.5 


(C) 


200 


16.5 


2.5 


(S) 


400 


12.2 


0 
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.1 : Simplified geological mnp of Ghorabi - Nasser area (simplified after El Aref and Lotfy, 1989) 

Clastic rocks of the Cenomanian Bahariya Fm.;2= karstified limestones of the Middel Eocene Naqb Fm. 
Iron ore deposits ; 4= Quaternary sediments; 5= Faults; 6= Drainage lines 
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Fig 2 . Zone Divisions 
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Fig 7. Directional Variogram for the North Zone (N). 
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Fig 11. directional Variogram Models for the Centeral Zone (C). 
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Fig 12. Directional Variogram Models for the North Zone (N). 
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Fig 13. Directional Variogram Models for the East Zone {E). 
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Fig. 14 Zonal Variogram Models. 
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